Ultrafast charge recombination in photoexcited Mott-Hubbard insulator 
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We present a calculation of the recombination rate of the excited holon-doublon pairs based on the two- 
dimensional model relevant for undoped cuprates which shows that fast processes, observed in pump-probe 
experiments on Mott-Hubbard insulators in picosecond range, can be explained even quantitatively with the 
multi-magnon emission. The precondition is the existence of the Mott-Hubbard bound exciton of the s-type. We 
find that its decay is exponentially dependent on the Mott-Hubbard gap and on the magnon energy, with a small 
prefactor which can be traced back to strong correlations and consequently large exciton-magnon coupling. 
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The advance in the ultrafast absorption spectroscopy allow 
for a direct access to excited states in a variety of materials 
as well as detailed studies of relaxation and thermalization 
processes. Correlated systems and Mott-Hubbard (MH) in- 
sulators are in this connection particularly of interest since 
here phenomena differ qualitatively from those well under- 
stood in band insulators and semiconductors. A prototype 
MH system are undoped cuprates as YBa2Cu3C<6, La2Cu04 
and Nd 2 CuC»4 which have been investigated with femtosec- 
ond pump-probe spectroscopy with an universal obser- 
vation of ultrafast relaxation and recombination. The photoin- 
duced carriers in undoped cuprates excited across the MH gap, 
which is in these systems of charge-transfer type, recombine 
in a picosecond range being a scale far below the usual physics 
of clean semiconductors with similar band gaps fl4£] . The ori- 
gin of such ultrafast processes has been already qualitatively 
attributed to strong-correlation effects, in particular to larger 
energy of magnon excitations with the exchange energy J and 
strong intrinsic coupling of charge carriers with spin fluctua- 
tions within the MH insulator yD, still no model-based theory 
was presented so far. Other MH insulators probed recently by 
ultrafast spectroscopy |0, Q] also exhibit fast-relaxation phe- 
nomena. Analogous are findings in the experiments with ul- 
tracold fermions in optical lattices where the decay of double 
occupancies (H 0] are related to the correlated nature of the 
background state. 

Due to strong Coulomb repulsion (strong correlations) ex- 
cited quasiparticles within the MH insulators can be con- 
sidered as empty sites (holons) and doubly occupied sites 
(doublons), respectively. In the following we consider as 
a prototype model the single-band Hubbard model on two- 
dimensional (2D) square lattice in the regime of large U ^> 
t, whereby the generalization to more appropriate charge- 
transfer model for cuprates is quite straightforward and the 
comparison with the parameters for actual materials can be 
also performed. 

When interpreting the pump-probe experiments it is as- 
sumed that initially photoinduced mobile charges - holons and 
doublons - after very fast transient in the femtosecond range 
form a MH exciton, i.e., a bound holon-doublon (HD) pair. 
The existence and stability of the MH exciton has been previ- 
ously addressed in connection with cuprates both in model 



studies II10I-I12I1 as well in the interpretation of optical re- 
sponse II7I J13I1 and of large Raman shift fl4ll with some con- 
flicting conclusions. As a result of our study the disagreement 
can be reconciled with the observation that the ground state 
(g.s.) of bound HD pair is of the s-type, thus not observable 
in an optical absorption. Just such a s-type pair state appears 
to be the precondition for a fast non radiative recombination 
process. The evidence for the intermediate HD bound state 
can be pump-intensity independent recombination rates, ap- 
parently observed in experiments Q3|] . Nevertheless, the for- 
mation of an MH exciton is a nontrivial consequence of strong 
correlations and large charge-magnon coupling in 2D antifer- 
romagnetic (AFM) background. The recombination of an MH 
exciton, i.e., a bound pair of a doublon in the upper Hubbard 
band (UHB) and a holon in a lower Hubbard band (LHB), 
across the MH gap A via the the multi-magnon emission is 
schematically presented in Fig. Q] Using the golden-rule ap- 
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Figure 1: (Color online) The process of the Mott-Hubbard exciton 
recombination via the multi-magnon emission. 

proach to the recombination rate T = 1/r of the exciton (nu- 
mericallly evaluated on small systems) we show that it is very 
large in comparison with clean semiconductors. In particular 
our results are consistent with the exponential dependence, 



T oc exp(— oA/J), 



(1) 



so that enhanced T has origin in large exchange scale J, but 
even more in small a ~ 0.5 which can be again traced back to 
strong correlations. In spite of large A recombination times 
t are in the picosecond range, consistent with pump-probe 
experiments in undoped cuprates 
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As the starting model for the study of photoexcited states in 
the MH insulator we use the single-band Hubbard model, 



U J2 



(2) 



where sums (ij) run over nearest-neighbor (n.n.) pairs of 
sites on a square 2D lattice. More realistic model for cuprates 
beyond the Hubbard model we comment furtheron. To re- 
duce the complexity and as well to get correlation effects 
more transparent in the regime U ^> t relevant for the MH 
insulator we perform the standard canonical transformation 
1 1511 . Within the lowest order in t 2 jU and neglecting the 3- 
site terms one gets a generalized t-J model for an insulator, 
including besides a spin-exchange term also both holon and 
doublon operators, 

H tJ = t^2{h] s h ls -d] s d ls +h.c.) + U^n di + 



(3) 



(ij) 



where h\ s — c. is (l — n is ) and d\ s = c\ s n is are projected 
fermionic operators, excluding double and unoccupied sites, 
respectively. Note that h\ s , d\ s up to a sign represent creation 
holon and doublon operators on a site previously occupied by 
spin s while n di = J2 S d\s d ™ and n h% = J2 S K s h ia are the 
local number of doublons and holons, respectively. It is ev- 
ident that H t j conserves the number of holons Nh and dou- 
blons Nd- For the study of the HD recombination it is essential 
to include the 3-site term, which emerges after the canonical 
transformation [15] also within the order t 2 /U. It has been 
so far mostly derived and considered for the hole-doped refer- 
ence insulator JlH 17 1 and can be for the HD case (neglecting 
terms which do not represent HD recombination) written in a 
spin-invariant form, 



H,=t 



hi. q d. 



lS > ■ Sj + h.c], 



(4) 



where i ^ k are n.n. sites to j and S the Pauli matrix. 
Within the lowest order of the t/U expansion one gets t$ = 
2t 2 /U = J/2. Evidently, H3 creates and annihilates HD pairs 
which are on next-nearest-neighbor (n.n.n.) sites. It should 
be stressed that within a canonically transformed Hamiltonian 
there is no recombination term being on n.n. sites or being 
linear in t. We use furtheron the hopping as the unit of energy, 
t = 1. 

We proceed to the evaluation of the HD recombination pro- 
cess in two steps. To describe the relaxation of the photoex- 
cited insulator through the intermediate bottleneck stage of an 
MH exciton we first establish the stability of the bound HD 
pair within the H t j, Eq. (f3]), neglecting H3. The problem is 
analogous to extensively studied case of the binding of two 
holes Nh = 2 within the t-J model 11811 . with an important 
distinction that holon and doublon are distinguishable parti- 
cles and can form also s-type (Ai symmetry) bound pair. In 



the Nh = 2 problem the g.s. has the d-type bound state above 
the threshold J > J c - 0.3 jfTsh . 

We study eigenstates of a single HD pair, i.e., Nh = Nd = 
1 by performing the exact diagonalization of H t ,j on small 
square and rectangular lattices with periodic boundary con- 
ditions using the Lanczos technique on N < 26 sites. As 
a criterion for binding we use the pair binding energy eb — 
E$ where E^ d , Eq, Eq, E$ refer to the 



jphd 
^0 



E$ 



HD pair, single hole, single doublon and undoped AFM g. s., 
respectively. Within H t j, Eq. (0, we also take into account 
Eq = Eq + U, while for 6q = Eq — Eq we use known ac- 
curate fits lH. Results for e 6 (J) obtained for N = 18, 20, 26 
are presented in Fig.|2]for the lowest s-type state. Values for e& 
appear to be quite size independent and confirm previous find- 
ings IU2I1 of the stability < of the MH exciton. The state 
is, however, of an even symmetry and not observable in an 
optical transition from the insulator AFM state. On the other 
hand, the p-type state relevant for optical absorption reveals 
e b ^ 0, hence does not appear to be a bound one. 
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Figure 2: (Color online) Holon-doublon pair binding energy e& vs. 
J. Plotted are results for the s-type g.s. for iV = 18, 20, 26 lattices. 



Another test for the HD binding are g.s. density correlations 
Dj = (tpQ \rihjndo\ipo ) (the position of the doublon chosen 
as the origin). In Fig. [3] Dj obtained (and symmetrized to 
recover the reflection symmetry) on N = 26 for J = 0.4 are 
consistent with the binding where the largest probability is for 
the HD pair being on a distance do = V^. On contrary, the 
unbound states (as the p-type lowest state) are characterized 
by largest do ~ V^V) for a given lattice. 

The next stage is to evaluate the recombination rate T = 
t -1 of the MH exciton using the Fermi golden rule with H3, 
Eq. (HJi, serving as the perturbation, 

T = 2^J2\^° m \ H ^o d )\ 2 KE° m - E% d ), (5) 

m 

and transitions are into magnon states |Vw within the un- 
doped AFM. We note that Eq. (0 can be represented as a re- 
solvent r = -2 ImCV = E' ld - £g), 



1 



H 3 \^ d ), (6) 



uj+ + E Q ' - Hj 
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0.4, as calculated for N 



with lu + = uj + iS. In the calculation only the AFM part 
Hj of the H t j, Eq. OJ, is relevant. Eq. (O is convenient 
for the evaluation within the Lanczos procedure, as usual for 
g.s. correlation functions it 1 8L I2QTI . Starting with the perturbed 
\4>q) = HalipQ ) the Lanczos procedure generates a tridiago- 
nal matrix and C{ui) can be expressed in a continued-fraction 
form. While 8 should be infinitesimally small thermodynamic 
limit N — > oo can be efficiently simulated by adopting finite 
S <C J. Moreover, to avoid the influence of such a Lorentzian 
smoothing (note that final r <C 1) we employ rather a Gaus- 
sian smoothing of poles produced after a finite number of 
Lanczos steps, M ~ 150 . 

The convergence of results with the system size is presented 
in Fig. [4] where T is evaluated for fixed J — OA and different 
lattices with N = 20, 24, 26 (N — 24 = 4 x 6 sites is a rect- 
angular lattice) sites. Results are shown versus the effective 
MH gap A = Eft* - between the bound HD-pair state 
and the AFM g.s. (see Fig. [T]). Taking A as a variable cor- 
responds to assume U in Eq. (f3]i as an independent parameter 
which could deviate from the relation U — At 2 / J obtained 
within the U t limit. We note from Fig. |4]that even small 
5 = 0.07 is enough to obtain continuous T(A) for the largest 
system N = 26. Clearly, the system size N should not be too 
small in order to accommodate enough magnons and to get 
reliable results in the relevant tails T <C 1. 




Figure 4: (Color online) Exciton recombination rate V vs. gap A 
for J = 0.4 and different planar systems with N = 20, 24, 26 sites. 
Results are smoothed with S — 0.07. 



As the central result we present in Fig. |5^,b calculated T 
vs. A obtained on the largest system N — 26 and different 
J = 0.3,0.4,0.6. The dependence on A/J in the relevant 
regime T(A ~ U) is close to the exponential, Eq. (HJ, as pre- 
sented in Fig.|5^. Effective a in this case is within the window 
0.3 < a < 0.7. Even better fit can be reached by following 
the perturbation-theory estimate for the probability of the gen- 
eration of rip bosons OSl IE IllD where in our case no = A/J, 



Figure 3: (Color online) Density correlations Dj of the hole with p q, 
respect to the doublon position for J 
26 lattice. 
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Results presented in Fig. |5]3 for T vs. (A/J)ln(A/te) reveal 
even better agreement, with quite universal cto ~ 0.8 (instead 
of ctQ = 2 following from straight derivation of Eq. (O ). 




6 10 
A/J ln{A/te) 
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Figure 5: (Color online) Exciton recombination rate Y vs. a) A/J 
and b) A/Jln(A/ie) for different J = 0.3,0.4,0.6 as calculated 
for N = 26 sites. 

Exponential dependence of V on J, Eq. ([TJ, is plausible 
since the recombination involves the multi-magnon emission. 
Taken at the face value within the reduced model, Eq. (0, 
with A = U - 2|ejj| - |e&| and U = At 2 / J the variation 
of r with J is even enhanced. Nevertheless, even for quite 
realistic value J = 0.4 we get in this simplified case A/J ~ 
17 and T ~ 10 -6 which is not unreasonably small taking 
into account large number of involved magnons. On the other 
hand, any decrease of the 'effective' gap A leads to a large 
enhancement of the rate V. 

Before the application of our theory and results to actual 
MH-insulator materials and experiments we offer some expla- 
nation of substantial T even for cases where evidently a large 
number of emitted magnons no ~ A/J^$> 1 is required for 
the HD-pair recombination. It is already evident from numer- 
ous numerical 111 811 and analytical [22] studies of the single 
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hole Nh — 1 within the t-J model that a holon (or doublon) 
perturbs strongly the AFM background in the strong correla- 
tion regime J < t. Even more this is the case for two (bound) 
holes Nh — 2 ll 811 and for the present example of a bound HD 
pair, Nh = Nd = 1. In Fig. |6]we present results for exchange- 
energy deviations relative to the reference AFM g.s. Results 
are for J = 0.4 and for the most probable HD configuration 
with the pair distance do = y/2. It is quite evident that the 
recombination of a HD pair via H3 (which just requires the 
holon and doublon being on n.n.n. sites) already leads to the 
generation of a large number of magnons no 1. 
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Figure 6: (Color online) Bond perturbations within the holon- 
doublon pair g.s., presented for the most probable pair configuration. 
Results are for J = 0.4 and deviations of bond energies relative to 
AFM g.s. are presented. 



Related but not equivalent argument for the appreciable V 
for cases A J is small parameter ao < 1 in the Ansatz 
and even more a ~ 0.5 (within the relevant regime of J/t) 
in Eq. 0. This can be considered as another manifestation of 
the strong coupling between the charge and spin excitations 
within the MH insulator as well as in strongly correlated sys- 
tems in general. 

Application to cuprates: Presented analysis and results 
can be generalized in quite straightforward way to insula- 
tors, where the MH gap represents charge-transfer excita- 
tions. In particular this includes insulating and superconduct- 
ing cuprates, where more complete (three-band) model in- 
cludes hybridized Cu d orbitals and O p orbitals. Although 
this goes beyond the prototype single-band Hubbard model 
in our study, through the Zhang-Rice mechanism [23] quasi- 
particles can be well described within the t-J model both for 
doped holes and doped electrons J24ll . Parameters are quite 
well established as t ~ 0.35eV and J/t ~ 0.4 [24]. The 
three site term H3, Eq. (0J, as well emerges from the three- 
band model lfl7ll . although £3 can be quantitatively different 
but also less important entering V in Eq. (0 only as a prefac- 
tor. So the essential generalization emerges from fact that the 
MH gap A and the 'effective bare' U in the model, Eq. 0), 
are not directly related to J and can be considered to large 
extent as independent parameters, taken from experiment or 
from more complete model. Optical MH gaps in undoped 
cuprates are experimentally known Ao ~ 1.5eV, still differing 
between members of cuprate family. In order to get defined 



A in Eq. and in Fig. we have to take into account that 
optical transitions are to HD p-type unbound states, so that 
A = A - je 6 |. 

We first consider Nd2Cu04 with the optical gap Ao = 
1.6 eV H and J = 0.155 eV, so that we get A = 4.1 t and 
r ~ 2.2 • 10 -2 in units of r^ 1 where r = h/t ~ 2 fs. This 
yields final relaxation time r = ~ 0.09 ps which is sur- 
prisingly close to experimentally established r ~ 0.2 ps J^t]. 
Analogous, for La2CuC>4 with Ao = 2 eV and J = 0.133 eV 
we get A = 5.3 t and T ~ 1.3 • 10 -4 which yields t ~ 15 ps. 
For the latter compound r has not been well determined [3] 
but is considerably longer than in Nd2Cu04 consistent with 
our result. 

Conclusions and discussion: The most important finding of 
this study is that multi-magnon emission can be quite an ef- 
ficient mechanism for the nonradiative recombination of pho- 
toinduced charges in a MH insulator. This is in contrast to 
usual band insulators and semiconductors where the impor- 
tance of analogous multi-phonon processes has not been es- 
tablished theoretically or experimentally [4], while such pro- 
cesses could be relevant for some novel structures as carbon 
nanotubes lf2~ll l25ll . The distinction of MH insulators is not 
just in larger boson scale J compared to phonon energies but 
primarily in the strong coupling between charge quasiparticles 
(holes and doublons) and spin excitations. The main manifes- 
tation is in the established prefactor ao < 1 in Eq. (0 as well 
as in a ~ 0.5 in the exponential Eq. (0). It should be also 
reminded that in spite of strong correlations doped holes and 
doublons are quite mobile with modestly enhanced effective 
mass. 

When interpreting recent pump-probe experiments on MH 
insulators, it should be reminded that we consider here only 
the bottleneck process of a MH-exciton recombination. The 
initial relaxation of the pump-induced transient metallic-like 
collective state of holons and doublons [3] is expected to be 
much faster and has been already investigated although not 
fully settled in different model studies [26-30]. 

It should be noted that a similar mechanism of the HD-pair 
recombination has been recently considered as a candidate for 
the relaxation in fermionic cold-atom systems JUlia] although 
it appears to be subdominant process to the kinetic-energy 
assisted decay which requires a finite density (metallic-like 
state) of nonequilibrium holons and doublons. The latter 
could emerge at effective T > or in systems far from equi- 
librium. On contrary, for 2D cuprates our study reveals that 
the multi-magnon mechanism of charge recombination and 
thermalization is very efficient and leads to ultrafast relaxation 
in the picosecond range as experimentally observed. 

The authors acknowledge fruitful discussions with T. To- 
hyama. This work has been supported by the Program Pl- 
0044 and the project J 1-4244 of the Slovenian Research 
Agency (ARRS). 
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